UNEXPECTED LOCAL EXTREMA 
FOR THE SENDOV CONJECTURE 



MICHAEL J. MILLER 



Abstract. Let S(n) be the set of all polynomials of degree n with all roots 
in the unit disk, and define d(P) to be the maximum of the distances from 
each of the roots of a polynomial P to that root's nearest critical point. In 
this notation, Sendov's conjecture asserts that d(P) < 1 for every P £ S(n). 

Define P £ S(n) to be locally extremal if d(P) > d(Q) for all nearby 
Q S S(n), and note that maximizing d(P) over all locally extremal polynomi- 
als P would settle the Sendov conjecture. 

Prior to now, the only polynomials known to be locally extremal were of 
the form P(z) = c(z n + a) with \a\ = 1. In this paper, we determine sufficient 
conditions for real polynomials of the form 



In 1958, Sendov conjectured that if a polynomial (with complex coefficients) has 
all its roots in the unit disk, then within one unit of each of its roots lies a root 
of its derivative. This conjecture has yet to be settled, although it has been the 
subject of more than 80 papers over the intervening years, and so has been verified 
for many special cases. These have been documented by Sendov [8], Schmeisser [BJ, 
Sheil-Small [9, Chapter 6] and Rahman and Schmeisser [5] Section 7.3]. 

Let n > 2 be an integer and let S(n) be the set of all complex polynomials of 
degree n with all roots in the unit disk. For a polynomial P with roots z\ , . . . , z n 
and critical points (f, . . . , Cn—i, define 



If P G S(n), then the Gauss-Lucas Theorem [5J Theorem 2.1.1] implies that each 
d(P) < 2, and Sendov's conjecture asserts that each d(P) < 1. 

We will say that a polynomial P is expected if it is of the form P(z) = c(z n + a) 
with |a| = 1. In 1972, Phelps and Rodriguez defined a polynomial P £ S(n) to be 
extremal if d(P) = sup{d(Q) : Q £ S(n)}, and conjectured [H after Theorem 5] 
that extremal polynomials are all expected. Since any expected polynomial P 
has d{P) — 1, this conjecture implies Sendov's conjecture. The Phelps- Rodriguez 
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conjecture has also been verified for a number of special cases, as documented by 
Rahman and Schmeisser [5j Section 7.3]. 

Define an e-neighborhood of P € S(n) to be the set of all polynomials Q £ S(n) 
whose roots are within e of the roots of P (in the sense that the roots of Q can 
be paired with the roots of P so that in each pair, the difference of the roots has 
a modulus less than e). Define a polynomial P s S(n) to be locally extremal if 
d(P) > d(Q) for all Q in some e-neighborhood of P, and note that maximizing 
d(P) over all locally extremal polynomials P would settle the Sendov conjecture. 

The expected polynomials are all locally extremal, as was demonstrated by 
Vajaitu and Zaharescu [10] and by Miller %, Theorem 3]. Given this, it is tempting 
to approach Scndov's conjecture by trying to show that all locally extremal poly- 
nomials must be expected. Indeed, Schmieder has made several such attempts [7], 
although Borcea has shown [21 section 1] that each contains a technical flaw. 

In this paper, we prove 

Theorem 1. For each n G {8, 9, 12, 13, 14, 15, 19, 20, 26}, there are locally extremal 
polynomials P £ S(n) of the form 



Note that polynomials of this form would not be expected, as each has a root ft that 
is not on the unit circle. Thus Theorem [1] implies that there are locally extremal 
polynomials that are not expected. 

Theorem[T]has several important consequences. First, it shows that the technical 
flaws in Schmieder's approach cannot be patched. Second, the structure of these 
polynomials (in particular, the multiple critical point) identifies circumstances that 
potential variational proofs of Sendov's conjecture will need to address. Finally, 
finding these polynomials is a significant new step in identifying all local extrema. 

In section [5J we list eight properties that are sufficient for a polynomial to be 
locally extremal, and in sections [3] and [4] we verify that these properties suffice. In 
section[5l we describe how to construct polynomials that satisfy all these properties, 
and in section [6] we list the resulting polynomials, thereby verifying Theorem [T] 



In section 4 we prove that for a real polynomial P of degree n > 5 to be locally 
extremal, it suffices for it to satisfy 8 properties, beginning with the following. 

A: All roots of P lie in the closed unit disk. 

B: All roots of P that are on the unit circle are simple. 

C: P has a root at ft, with < (3 <1. 

D: All critical points of P lie on a circle of positive radius centered at ft. 

E: P has a real critical point a < [3 of order n — 3. 

Let Z\, . . . , z n be the roots of P, numbered so that Z\, . . . , z m are on the unit 
circle. Let (i, ■ ■ ■ > Cn— l be the critical points of P, numbered (as allowed by prop- 
erty [E| so that Q = a for j > 3. 

Note that property [Dl implies that P has a simple root at [3. Our next property 

is 

F: For each Zi ^ ft, we have minix^xn-i \zt — Cj\ < mini<j<n-i \ft ~ Cj\ < 1- 




(w-a) n 3 (w 2 + bw + c) dw with < ft < 1. 



2. Properties 
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To examine the effects of changing and the Q by small amounts, denoted by 
A/3 (which we will require to be real) and A(j, we will use the following notation. 

"AC fc -A/3l 



(2.1) 



Ek 
F k 

E 3 



-3? 



for fc = 1 and k = 2, 
»ES AO] - (n - 3)A/3 



for k = 1 and fc 



/3-a 

e;^(3[ao]) 2 



and 



2(/3-a) 2 ' 

Since (up to a constant multiple) P'(z) = YVj=i( z ~ 0) an d by property ICl we 
have P(z) = f* P'(w) dw, then the roots of P are functions of /3 and the Q. By 
property[Bj the roots of P that are on the unit circle (being simple) are differentiable 
functions of /3 and the Q . 

Recall that z\,...,z m are on the unit circle. For i = 1, . . . , m dehnc 



E, 



1 dzj 
Zi d/3 



A/3- 



i=i 



1 dzj 



i+3 





f z > P'(w)dw~ 


^ZiP'iZi) J 


13 (w- a) 2 _ 



AC, 



n-l 

E 

3=3 



(9[A0 



1 ggj 

2* d( 3 



i=3 



and 



Note that for a fixed polynomial P, the m + 3 expressions Ek are all linear 
in the 7 real "variables" A/3, K[AO], 3[A&], 3?[AC 2 ], 3f[A<a], ^E"=3 A Q] and 
^[^""g 1 AO], and that the m + 3 expressions are all constant multiples of the 
real "variable" S™ = 3 (S[A0]) 2 - Our final two properties are 

G: There are constants Ck > (depending on P, but independent of the 8 "vari- 
ables") so that the sums ' c k E k = and YJk=i c kF k = E^^AO]) 2 . 



H: The coefficient matrix of the system {Ek = 
"variables" is of rank 7. 



k = 1, . 



, , m 



3} in our 7 



3. Preliminary calculations 

We will show that a real polynomial P of degree n > 5 satisfying properties EHH 
is locally extremal, as follows. 

Define r to be the radius of the circle in property O and recall that (up to a 
constant multiple) P(z) = Jp YVj=i(' w ~ ) dw. From our properties \K\ and [Dl we 
know that P has all its roots in the closed unit disk and all its critical points on a 
circle of radius r > centered at /3. Define the n-tuple (A/3, AO, ■ ■ ■ , AOi-i) € C" 
to be an improvement of P if A/3 is real and if the polynomial 

, z n— 1 

(3.1) / n>-(& + A C;)]^ 

■//3+A/3 J = 1 

has all its roots in the closed unit disk and all its critical points strictly outside the 
circle of radius r centered at (3 + A/3. By property IP1 we know that r < 1, so there 
is at least one improvement, namely (1 — /3, — 0> • • • , — Cn-i)- 
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For an improvement / = (A/3, A£i, . . . , A£ n -i), define 

/ n-l 

||/||= |A/3| 2 + ^|A0| 2 



Note that ||/|| > 0, for if ||/|| = 0, then the critical points of (|3.1j) would be on 
(and thus not strictly outside) the circle of radius r centered at ft + A/3. 

Take any e with < e < ft/2 and recall our definition of e-neighborhood from 
section[TJ The critical points of P are continuous functions of the roots of P, so there 
is a S > such that for all polynomials Q in a ^-neighborhood of P, the roots and 
critical points of Q are within e of the roots and critical points of P. By PropertvlF] 
we know that d(P) = r, so if P is not locally extremal, then there is a polynomial 
P G S(n) in the 5- neighbor hood of P with d(P) > r. Let ft be the root of P paired 
with ft (so |/3— ft\ < e), and for j — 1, . . . , n— 1, let Q be the critical point of P paired 
withCj (so \C 3 -C\ < e)- Note that \ft\ > \ft\-\ft-ft\ > ft/2 > 0, define u = |/3|//3and 
note that \u — 1 = ||/3| — ft\/\ft\ < 2e/(ft/2). Since \u\ = 1, then the transformation 
z — > uz is a rotation about the origin, so / = (uft — /3, uCi — £i, . . . , uCn-i — Cn-i) 
is an improvement of P. Note that each — Cil < \ u \\(j ~ 01 + 10 1 ( u — 1) — C £ 
for some constant C (and similarly for \uft — /3|), so ||J|| can be made arbitrarily 
small. 

We have just seen that if P is not locally extremal then there are improvements 
I of P with ||/|| arbitrarily small. This means that we can show that P is locally 
extremal by proving 

Theorem 2. If P is a real polynomial of degree n > 5 that satisfies properties \A\U^ 
then there is a constant A > so that every improvement I of P has \\I\ \ > A. 

We begin with some preliminary calculations. For t > 0, we'll say that a (real 
or complex) quantity is 0(t) if its modulus is bounded by a constant multiple of t. 
Given this, we have 

Lemma 3. If z = 0(t) then 

(1) \l + z\ = 1 + 5R[z] + 0(t 2 ), and 

(2) |1 + z\ = 1 + Sl[z] + (l/2)(3[z]) 2 + 0(t 3 ). 

Proof. Note that $l[z] = 0(t) and \z\ 2 = 0(t 2 ). The results follow from the equality 
|l + z| = (1 + 2^^ + lzl 2 ) 1 / 2 , the Taylor series (1 + s) 1 / 2 = 1 + s/2 - s 2 /8 + 0{s 3 ) 
and the substitution s = 2^R[z] + \z\ 2 . □ 

We will be examining the relationships between the roots and critical points 
of P. We calculate the partial derivatives of these relationships with 
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Lemma 4. Write P'(z) = Y[j=i( z ~ 0) an d suppose that z% is a simple root of 
P(z) = f*P'(w)dw. Then 

(1) Zi is an analytic function of P, Ci, ■ ■ ■ , Cn-lj 
(2) 



dzi 


P'(P) 




d/3 ~ 






dzi 


1 r Zi P'(w)dw 




dQ ~ 


P'(zi)J w-Q ' 




d 2 z t 


2 dm P"(zi) 


fdzA 2 


dQ 


Zl - Q dQ P'(z t ) 


\dQJ 


d 2 Zi 


1 dzt 1 


dz l 


dCjdCk 


z% - Cfe d(j Zi - 


Ci d( k 




P"{ Zi ) & Zl d Zl 


1 



(4) ^ = ^^7^--^Ft[7^) '""I 
(5) 



P'(w)dw 

for j f= k. 



P>( Zl )dQd( k P'{ Zl )J p {w~Q){w-Ck) 

Proof. We may assume that Zi ^ (3 (else the results would be trivially true.) Proofs 
of parts 1-3 can be found in [TJ Lemmas 2.1 and 2.3]. Part 4 can be established by 
writing P'(z) = (z — (j)Q(z) (with Q independent of Q) and calculating the second 
partial derivative (with respect to Q) of 



= / wQ(w) dw — Q / Q{w) dw. 

Part 5 can likewise be established by writing P'(z) = (z — Cj)i z — Ck)Q(z) and 
calculating the mixed second partial derivative. □ 

4. Proof of Theorem 

Recall that P is a real polynomial of degree n > 5 that satisfies properties [XI IHl 
We will prove Theorem [2] in two stages, first by estimating the roots and critical 
points of P with linear approximations, and then by improving our estimates with 
quadratic approximations. 

We begin by examining the conclusions that can be drawn from linear approxi- 
mations to the roots and critical points of P with 

Proposition 5. If I — (A/3, A£i, . . . , ACn-i) is an improvement of P and At — 
\\I\\, then A/3 = O(At) 2 , AQ = O(At) 2 for j < 2, R[AQ] = 0{At) 2 for j>3 and 
E"=3 A0 = O(At) 2 . 

Proof. From the definition of At we know that A/3 = O(At) and that each AQ = 
O(At). Since I is an improvement, then each |(Cj — 0) + {AQ — A/3)| > r, so 
1 + (ACj — A/3)/(Cj — /3)| > r/|Cj — /3| = 1, so using part 1 of Lemma [3] gives us 
that each 



AC, - A/3 



b-0 



< O(At) 2 



(4.1) -SR 

Recalling the expressions E k defined in 12.11 this gives us two inequalities E\ < 
O(At) 2 and E 2 < O(At) 2 . Adding igU} for j = 3, ... ,n - 1 and recalling that 
Cj = a is real for j > 3 gives us a third inequality E3 < O(At) 



2 
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By part 1 of Lemma [4] the roots of our improved polynomial (13. 1| that originate 
from the simple roots z% of P are analytic functions of j3, Ci, . . . , ( n -i, so each such 
root is of the form Zi + Azi, with 

Note that each Az, — O(At). Since / is an improvement, then each \zi + Azi\ < 1. 
If \zi\ — 1, then |1 + Azi/zi\ < l/|z;| = 1 so using part 1 of Lemma [3] gives us 
^[Azi/zi] < O(At) 2 and thus we have inequalities E i+3 < O(At) 2 for 1 = 1, ... , m. 

By property [Gl there are positive constants Ck so that Y^k=i c kEk — 0. Since 
each Ei < O(At) 2 , then each 

m+3 

E t . = (-I/q) C * E * > O(At) 2 

so each Ei = O(At) 2 . 

Thus we consider the system {E^ — O(At) 2 : k = 1, . . . , to + 3}. By property [Hi 
the coefheient matrix of this system is of rank 7, and so solving this system shows 
that the values of our 7 "variables" are all O(At) 2 . Thus we can conclude that 
A/3 = 0{At) 2 , that AQj = 0{At) 2 for j < 2, and that Y%=3 A = 0{At) 2 . 

Suppose that j > 3 and note that J2k=l ^[^Cfe] = O(At) 2 . Since by property|E] 
we know that Q — (3 — a — (3 < 0, and since A/3 = O(At) 2 , then (|4.1|1 implies that 
each sft[A£j] < O(Ai) 2 , so each 

5R[A0] - - ^[ A Cfe] + O(Ai) 2 > O(Ai) 2 

j#fe=3 

and so each 9fJ[A(*,-] = O(At) 2 . This finishes the proof of Proposition [5l □ 

At this point, for j > 3 we know only that each AQ = O(At). We can improve 
this estimate by looking at quadratic approximations to the roots and critical points 
of P with 

Proposition 6. If I — (A/3, A£i, . . . , ACn-i) is an improvement of P and At — 
then each AQ = O(At) 3 / 2 . 

Proof. Note that the hypotheses of Proposition O are satisfied, so we may use all 
of its conclusions. In particular, we know that AQ = O(At) 2 for j < 2 and that 
5R[A(*j] = O(At) 2 for j > 3, so to verify Proposition [5] we need only show that 
3[A(*,-] = O(At) 3 / 2 for j > 3. We will do this by repeating the calculations of 
Proposition but working now to O(At) 3 . 

Since / is an improvement, then each |1 + (AQ — A(3)/((j — (3)\ > 1, so using 
part 2 of Lemma [3] gives us 

2 



(4.2) - 5ft 



AC, - A/3 



AC, - A/3 

0-/3 



< O(At) 



Define inequalities 1 and 2 by evaluating (|4.2[) for j = 1 and j = 2 respectively. By 
Proposition [5l for j < 2 we have A^ — A/3 = O(At) 2 , so inequalities 1 and 2 are 
E k + F k < 0{Atf for k = 1, 2. 
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Define inequality 3 to be the sum of (|4.2| for j = 3, . . . , n — 1. Now each C,j—(3 = 
— ((3 — a) is real for j > 3, and each ^[ACj — A/3] = Q^A^], so inequality 3 can be 
written as E 3 + F 3 < O(At) 3 . 

Each root of our improved polynomial is of the form Zi + Az{. A quadratic 
approximation to Azi includes terms of the form A(3AQ = O(At) 3 and (for j < 2 
or k < 2) AQA(k — O(At) 3 , which are absorbed into the O(At) 3 when we write 

n ^ Q ^) n — 1 n — 1 n — 1 ~2 

^ - %>* + E + §§ E ^ + 5 E E + 

Note that Proposition [5] implies that each Azi = O(At) 2 . Now (j = a for j > 3, so 
Lemma |4] shows that for j > 3 and fc > 3 we have 

92z * 92% = Jp'(z i )A (™-«) 2 J 



From Proposition [5] we know that 



if jjtk. 



j=3 fc=3 \j=3 / 

and that (AC 3 ) 2 = -(3[ACj]) 2 + 0(Ai) 3 for j > 3, so 

n-ln-1 a2 r n-ln-1 , „ 2 „ 2 „ \ 

E E ^ ao ^ -EE(^fe-^ a & A & + ^ 



j=3 fc=3 ,J ' 3=3 fc=3 

n-1 



1 /" Zi P'{w)dw \ . , . ., , , 
' ^(3[A0]) 2 + O(At) 3 



i=3 



Recall that each \zt + Azi\ < 1 and that each Azi = O(At) 2 . If Zi| = 1, then 
1 + Azi/zi\ < 1, so using part 1 of Lemma [3] gives us $l[Azi/zi] < 0(Ai) 4 and so 
E i+3 + F i+3 < O(At) 3 for i = 1, . . . , m. 

Thus we have Ek + Ff. < O(At) 3 for k — 1, . . . , m + 3. From property[G] we have 
positive constants Ck so that ^"^(^[ACj]) 2 = Y^k=i c k(Ek + Fk) < O(At) 3 and 
thus 3[ACj] = O(At) 3 / 2 for j > 3, which completes the proof of Proposition [6l □ 

Using our estimates from Propositions [3] and [SI we can now write the 

Proof of Theorem^ Let I = (A/3, AO, . . . , A£„_i) be any improvement of P, and 
let At = ||i||. From Propositions [5] and we know that A/3 = O(At) 2 and that 
each A(j = £>(At) 3 / 2 , so 



A* = ||/|| = |^|A/3| 2 + | AC, | 2 I = O(At) 3 / 2 
Thus there is a constant if > so that At < K(At) 3 / 2 , so ||J|| = At > 1/K 2 . □ 
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5. Calculations 

Recall that Theorem [2] implies that any real polynomial of degree n > 5 that 
satisfies Properties lATlHl is locally extremal. In this section, we show how to find lo- 
cally extremal polynomials by describing how to construct polynomials that satisfy 
our properties lATlHl Maple code for these calculations can be found in the source 



files for this paper, located at http://arxiv.org/abs/math/0505424 [v3]. 



Recall that P is to be a real polynomial of degree n > 5. By property lEl we 
must have P'(z) = (z — a) n ~ 3 (z 2 + bz + c) for some real numbers a, b and c. 
By property [Cj we know that P(z) — J g P'(w)dw, so the coefficients of P arc 
polynomials in {/3, a, b, c}. 

Applying property [D] to the critical points of P (the real number a and the 
complex roots of z 2 + bz + c), we get our first equation 1 + bf3 + c — ([3 — a) 2 . 

Note that property iGl implies that the rows of the linear system {Ep. = : k = 
1, . ..,rn + 3} are linearly dependent. Since property [H] states that the coefficient 
matrix of this system is of rank 7, this means that there must be at least 8 equations 
in the system, so m > 5 and thus we must look for polynomials with at least 5 roots 
on the unit circle. 

If n is odd, we will look for polynomials with three pairs of complex conjugate 
roots on the unit circle. Each such pair will be the roots of a quadratic of the form 
z 2 +diZ + l, so the remainders upon dividing P by each of these three quadratics are 
linear polynomials with both coefficients equal to 0. This generates an additional 
six equations in the seven variables {/?, a, 6, c, di, c?2, cfe}. Thus we have a nonlinear 
system of 7 equations in 7 unknowns. 

If n is even, we will look for polynomials with two pairs of complex conjugate 
roots on the unit circle and a root at —I. The conjugate roots generate an addi- 
tional 4 equations in the six variables {/3, a, 6, c, d±, c?2} (as above), and the equation 
P(— 1) = generates a sixth equation. Thus we have a nonlinear system of 6 equa- 
tions in 6 unknowns. 

Thus in cither case we get a nonlinear system of equations, with the same number 
of equations as unknowns, so we can try to solve this system. Note that there may be 
more than one solution, so we will need to choose the "correct" one. For each such 
solution, we verify properties lATfHl by numerically checking the following assertions. 
(Details of these computations may be found in the Maple code referenced above.) 

Property [Aj The maximum modulus of the roots of P is equal to I (to the 
accuracy calculated.) 

Property [Bj The minimum distance between any two roots of P is greater 
than 0.1. (This shows that all roots of P are simple.) 

Property O We have 0.7 < < 0.9. (We know that P(f3) = 0, since P(z) = 
Jo P'(w) dw by construction.) 

Property [DI The distances between (3 and the critical points of P are all equal 
(to the accuracy calculated), and this common distance is greater than 0.9. 

Property lEl The quantity (3~a > 0.9. (Note that the critical point a is real and 
of order n — 3 by construction.) 

Property [F] If we define 



R = max < min |zj — Cjl : z i 7^ P ( an( l r = mm \fl — 

I l<j<n— 1 J 1 

then r < 0.97 and r - R > 0.02. 
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Property [Gj The linear system given by the sums has a solution in which every 
c k > 0.3. 

Property|H) The seventh largest singular value of the coefficient matrix is greater 
than 0.04. 

Once we have verified properties [XriHl for a specific polynomial P, we know by 
Theorem [2] that P is locally extremal and we are done. 

6. Proof of Theorem Q] 

For the values of (3 and P'{z) given below, one can verify that the polynomials 
P(z) = fp P'(w)dw satisfy Properties IAhTII and thus by Theorem [2] are locally 
extremal. (Details of these computations may be found with the Maple code refer- 
enced above.) 

For n = 8, we take (3 = 0.7290857513 and 

P'{z) = {z + 0.2035409790) 5 (z 2 - 0.5410836525z + 0.7327229666). 
For n = 9, we take = 0.7145672829 and 

P'{z) = {z + 0.2157115753) 6 (z 2 - 0.8021671918z + 0.9280147829). 
For n = 12, we take (3 = 0.8403619619 and 

P'{z) = (z + 0.1155828545) 9 (z 2 - 0.4090272613z + 0.5513532168). 
For n = 13, we take (3 = 0.8275325585 and 

P'(z) = (z + 0.1246203379) 10 (z 2 - 0.5415308686^ + 0.6699194279). 
For n = 14, we take (3 = 0.8158105092 and 

P'(z) = (z + 0.1304708647) 11 (z 2 - 0.6885970233z + 0.7916663399). 
For n = 15, we take (3 = 0.7999767588 and 

P'(z) = (z + 0.1400336168) 12 (z 2 - 0.8389864647z + 0.9148263642). 
For n = 19, we take (3 = 0.8684432238 and 

P'{z) = {z + 0.0923361850) 16 {z 2 - 0.6503807257z + 0.7337221736). 
For n = 20, we take (3 = 0.8570396874 and 

P'O) = {z + 0.0982636528) 17 {z 2 - 0.7563752823z + 0.8263310816). 
For n = 26, we take f3 = 0.8817716692 and 

p'(z) = ( z + 0.0797127446) 23 (z 2 - 0.7969496845z + 0.8496586550). 
These verifications complete the proof of Theorem [T] 
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